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ABSTRACT 

We present combined VLA and GBT images of NH3 inversion transitions (1,1) and (2,2) toward 
0MC2 and 0MC3. We focus on the relatively quiescent Orion cores, which are away from the 
Trapezium cluster and have no sign of massive protostars nor evolved star formation, such as IRAS 
source, water maser, and methanol maser. The 5" angular resolution and 0.6 kms -1 velocity resolution 
of these data enable us to study the thermal and dynamic state of these cores at ~ 0.02 pc scales, 
comparable to or smaller than those of the current dust continuum surveys. We measure temperatures 
for a total of 30 cores, with average masses and radii of 11 M© and 0.039 pc, respectively. Compared to 
other Gould Belt dense cores, the Orion cores have an unusually high gravitational-to-kinetic energy 
ratio (virial mass ratio R V i r » 1), resembling results for other clouds forming high-mass stars. This 
results from Orion cores having velocity dispersions similar to those in, e.g., Perseus and Ophiuchus, 
but higher masses for given sizes. 12 out of 30 cores are associated with embedded YSOs, which are 
identified by Spitzer. These Prostellar cores have a mean mass of 16.5 M versus 7.3 M Q for that 
of the starless cores. 14 of starless massive Orion cores are supercritical (mass-to-critical-mass ratio 
R c > 1). These massive Orion starless cores will likely collapse or fragment relatively quickly and can 
thus be considered direct precursors to protostars. 

Subject headings: ISM:clouds - individual (Orion) - Stars:formation - instrumentation: interferome- 
ters 



1. INTRODUCTION 

Stars form in molecular clouds. Within molecular 
clouds, discrete structures with observable column den- 
sity contrast, particularly in high density tracers, with 
its surroundings are referred to as cores (e.g. Benson & 
Myers 1989; Ward-Thompson et al. 2007). The relatively 
high extinction and density of cores make them the likely 
site for the onset of collapse, which could lead to star for- 
mation. 

Recent core surveys (e.g. Ikeda et al. 2007; Konyves 
et al. 2010; Sadavoy et al. 2010) are facilitated by focal 
plane imaging arrays of growing sizes. Although many 
surveys have large sample size in the hundreds, the ma- 
jority of these studies are of dust continuum with one 
of the major focus being the core mass function. Direct 
measurement of cores' thermal and dynamic structures 
require spectra maps of high density tracers and prefer- 
ably cover the same spatial dynamic ranges as those of 
the dust emission. 

There is a clear dichotomy between low mass star for- 
mation and high mass star formation. Massive stars are 
formed exclusively in giant molecular clouds and with 
higher efficiency. Due to the large distances of most 
massive star forming regions, many massive cores are 
under-resolved, showing signs o f star formation well un- 
derway, such as H2O masers (jMookerjea et al. I l200l 
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and/or compact HII regions (IReid fc Wilson 1 12005 K At 
about 437 pc (Hirota et al. 2007; Menten et al. 2007), 
Orion molecular clouds are the closest massive star form- 
ing region with an OB cluster, thus particularly suited 
for studying the early stages of star formation in massive 
cores and/or under the influence of young massive stars. 
In a series of papers, we identified quiescent clouds/cores 
(containing no HII region, no IRAS point sources, and 
are at least 1 pc away from the OB cluster) in Orion 
with NH3 and N2H" 1 " with a beam size of about 1' (Li, 
Goldsmith & Menten 2003; hereafter paper I), presented 
high resolution 350 /im images with a beam size of 9" 
(Li et al. 2007; hereafter paper II), and revealed that 
two thirds of Orion cores have signatures of ongoing dy- 
namic evolution, both outflows and inflows (Velusamy et 
al. 2008 hereafter paper III). Based on dust mass and 
dust core size (paper II), the majority of the cores are 
seemingly supercritical, i.e., no adequate support from 
thermal pressure, turbulence, or static magnetic field. 
This is consistent with the majority of the cores being hy- 
drostatically unstable (paper III), which, however, could 
not be directly tested due to a lack of high resolution 
spectroscopic data. 

Spectroscopic survey of CO 1-0 (Williams, Plambeck 
& Heyer 2003) and CO 3-2 (Takahashi et al. 2008) reveal 
that part of the Orion molecular cloud, namely, OMC2 
and OMC3, contains many molecular outflows. OMC2 
and OMC3 have also been mapped in high density trac- 
ers. Using the Nobeyama 45m telescope, Tatematsu et 
al. (2008) identifed 34 cloud cores in N 2 H+ 1-0 with a 
beam width of around 18". Ikeda et al. (2007) studied 
dense gas in the same region using H 13 CO + also with 
Nobeyama. These studies find higher column density for 
respective tracers compared with low mass star forming 
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regions, such as Taurus. NH3 inversion transitions are 
particularly suited for studying dense cores due to their 
lack of depletion and their sensitivity to kinetic temper- 
ature (Ho & Townes 1983). Wiseman & Ho (1998) ob- 
tained NH 3 maps of a 8' x 6' region around Orion-KL 
using VLA. The gas morphology there is dominated by 
quasi-parallel filaments severely influenced by the energy 
output of young massive stars. 

In this letter, we present the combined VLA and GBT 
ammonia survey of the cores in OMC2 and OMC3. The 
spatial resolution of ~ 5" and the spectral resolution of 
0.6 km/s facilitate a detailed examination of the thermal 
and dynamic states of massive 'quisecent' Orion cores. 

2. OBSERVATIONS AND DATA REDUCTION 

The VLA observations of OMC2/OMC3 were carried 
out in the D configuration on July 29 and September 24 
of 2000. We used the correlator mode 4 to cover the NH3 
(J,K) = (1,1) and (2,2) inversion transitions simultane- 
ously. With a primary beam of about 2' at the observing 
frequencies, a total of 20 pointings were used to mosaic 
the OMC2/3 regions. The correlator was configured to 
a bandwidth of 3.13 MHz for each transition that was 
divided into 64 channels, providing a velocity coverage 
of 40 kms -1 in 0.6 kms -1 channel spacing. In both 
observations, we used quasar 3C286 as the flux calibra- 
tor, 3C84 or 3C273 as bandpass calibrators, respectively. 
The time dependent gains were monitored by observing 
0539-057 or 0605-085. Visibility data were calibrated in 
AIPS and exported to MIRIAD for further imaging. The 
lcr RMS in the channel maps, after combining the two 
observations is about 8 mJy per spectral channel. 

The observations of NH3 (1,1) and (2,2) transitions 
with the GBT were taken on April 6th and 7th of 2005 
in the OTF mode using the receiver Rcvrl8-26. The cor- 
relator setup had a 200 MHz bandwidth and 0.31 kms -1 
velocity resolution (8192 channels). The flat baseline and 
the wide bandwidth allowed us to use frequency switch 
with a throw of 12.5 MHz. The spectra were folded and 
the 3D data cube were resampled to a grid of 12" spacing 
and 0.3 kms -1 channel. The resulting GBT spectra has 
a typical 0.12 K RMS noise per channel. 

We combine the NH3 data from the VLA with those 
from the GBT to recover missing short spacing fluxes 
in the interferometric data. The combination was per- 
formed in the UV domain using MIRIAD, following the 
procedure outlined in Vogel et al. (1984) and Zhang et 
al. (2000). The NH3 emission from the combined image, 
when convolved to the 30" GBT beam, recovers more 
than 80% of the fluxes detected in the GBT data. The 
integrated intensity of the NH3 (1,1) line with previously 
identified continuum sources and a pair of typical spectra 
are presented in Fig. [T] 

3. DERIVATION OF KINETIC TEMPERATURE AND 
VELOCITY DISPERSION 

Largely following Ho & Townes (1983), Paper I de- 
scribed a recipe for deriving kinetic temperature from 
spectrally resolved, modestly blended NH3 lines (intrin- 
sic line width AV ~ 1.0 kms -1 ). In this paper, we 
report line width AU in FWHM. The one dimensional 
velocity dispersion a used in Eq. [8l Eq. [TT] and Eq. [12] 
is related to FWHM as AV = i/8Tn(2)<7 for a Gaus- 
sian. The key step is to obtain the optical depth of the 



(1,1) line from simultaneously fitting all hyperfine com- 
ponents. As can be seen in Fig. [1] the VLA band is not 
wide enough to fully cover the two outer groups (osg). 
The combined data set thus only have the main and in- 
ner satellite components of the (1,1) transition and have 
a velocity resolution of 0.6 kms -1 . Toward most posi- 
tions, the intrinsic line width is smaller or approaching 
channel width. For such under-resolved line, the opacity 
cannot be easily uniquely fitted. 

We developed a more straightforward recipe 
utilizing two ratios between integrated inten- 
sities, which are directly observable, 1Z 12 — 

[jT^ na+is9) (l ) l)du]/[jT A (2,2)du] and TZ sm 

[fT^ S9 \l,l)dv]/[jT^ n9) (l,l)du]. The rotational 
temperature can be derived as the following 

T R = 41.5 K/ In [1.06 x C(l, 1) x R 12 } , (1) 

where C(l, 1) is a correction factor numerically deter- 
mined as 

C(l,l) = 0.003 + 2.26 K sm + 0.00032 e 538KSm , (2) 

which is based on fitting to simulated NH3 (1,1) spectra 
generated with a grid of opacities and linewidth. The 
kinetic temperature is then (paper I) 

T k = 3.67 + 0.307 xT R + 0.0357 x T| . (3) 

The full recipe is given in Kauffmann et al. (2012) and 
is generally applicable to a wide range of conditions in 
opacities, channel width, and intrinsic linewidth. The 
uncertainty in derived kinetic temperature is then calcu- 
lated using a Monte Carlo approach (paper I). For pixels 
with (2,2) detection, the representative lcr uncertainty 
is about 1 to 2 K for emission peaks and between 2 K 
and 5 K for diffuse areas. Except for about 1.5% of scat- 
tered pixels, the derived temperatures are between 10 K 
and 30 K. When there is no detection of the (2,2) line, 
we derive an upper limit to kinetic temperatures (T£) 
assuming 3ct(2,2) intensity. Figure [2] shows the derived 
temperatures overlaid with the dust continuum and dust 
cores. The temperatures are generally lower toward the 
inside the dust cores, which has been demonstrated in 
paper I. A quantitive analysis of the thermal structure of 
the whole region is presented in Kauffmann et al. (2012). 

The low velocity resolution of 0.62 kms -1 complicates 
the calculation of velocity dispersions. We first iden- 
tify main group channels with signal-to-noise ratios > 3 
are identified. If more than one such channel is found, 
a dispersion Cdata is calculated, in which channels are 
weighted by their intensity; if a single channel is found, 
an upper dispersion limit of 0.26 kms -1 is adopted. The 
intrinsic dispersion of the (1,1) main group is then sub- 
tracted to obtain er obs = [o^ ata - (0.2 kms" 1 ) 2 ] 1 / 2 . The 
value of the correction factor is derived from numeri- 
cal simulations of lines with various optical depths. The 
pixel averaged velocity dispersion for dust cores are listed 
in Table 1. 

4. DYNAMIC STATE OF CORES 

For a full coverage of all OMC2 and OMC3 cores, we 
use the SCUBA survey of Orion by Nutter & Ward- 
Thompson (2007), who identified cores as ellipses in the 
signal to noise ratio map. We select SCUBA cores with 
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Fig. 1. — The integrated intensity of NH3 (1,1) overlaid with dust continuum peaks from Chini et al. (1997) and references therein. The 
white dotted line indicates the location where the gain drops to 80% of the peak value. The spectra shown on the right are for the FIR 4 
position, the brightest NH3 peak in OMC 2. Due to the limitation of the VLA bandwidth, only the main group (mg) and the inner satellite 
groups (isg) of the hyperfine components are fully covered. The outter satellite groups (osg) arc on the band edge. 

more than 50% of the pixels having a measured TJt based 
on NH3. These temperature measurements are weighted 
by dust column density and then averaged to be the core 
temperature T core . Nutter & Ward- Thompson (2007) 
calculated dust mass assuming a uniform dust temper- 
ature of 20 K. In our map, 74% of the pixels have de- 
rived temperature lower than 20 K. We therefore scale 
the core mass by the factor -Bs50,20 k/ -E>850,T core , where 
Bs50,t is the Planck function at 850 /im and temperature 
T. The resulting core mass is generally larger than what 
is reported in Nutter & Ward-Thompson (2007). The 
derived core temperatures and core velocity dispersion 
are reported in Table [T] 

We identify cores with and without embedded stars 
based on Spitzer source identifications provided by 
T. Megeath. A core is deemed protostellar if the distance 
between its peak position and a YSO are less than the 
semi-minor axis of the core ellipse (Fig. [2] and Table [J). 



4.1. 



Are Cores Gravitationally Bound? 

McKee & Zweibel (1992) derive the virial theorem in 
Eulerian form to express explicitly the contribution of 
turbulent pressure both inside and surrounding the cloud 



is determined by the ratio of the difference between sur- 
face turbulent energy £q w and internal turbulent energy 
£ w , and total surface kinetic energy £q and total internal 
kinetic energy £ = |Afer 2 . For a cloud with primarily 
turbulent motions, c w = 3/2. If steady motions, such as 
rotation, are significant, c w is close to unity. The veloc- 
ity resolution of our data limits our ability to probe the 
thermal linewidth. The existence of blending of hyper- 
fine components in some spectroscopically resolved pix- 
els (apparent AV > 0.6 kms" 1 ) suggests that the gas 
there is unlikely to be in a pure sonic state of either a 
H2 thermal linewidth of AV ~ 0.3 kms -1 or even a NH3 
thermal linewidth of ~ 0.08 kms" 1 . For a large fraction 
of the pixels, however, we only have an upper limit to 
the linewidth. A unity c w is assumed for simplicity for 
the virial mass calculation in this section. 



Dust 



Kinetic T 



-5:00:00 



-5:06:00 



Q -5:09:00 - 





5h35m36s 5h35m21.6s 5h35m7.2s 
RA(J2000.0) 



5h35m36s 5h35m21.6s 5h35m7.2s 
RA(J2000.0) 



Fig. 2. — Left: 850 fim dust images with fitted dust cores (Nut- 
ter &c Ward-Thompson 2007). The red stars are embedded YSOs 
based Spitzer IRAC and MIPS surveys (T. Megeath, private com- 
munication) Right: Derived kinetic temperatures (no upper limits) 
overlaid with dust cores. Only cores with sufficient detection of 
temperature or upper limits (see discussion in the text) are plotted 
in this panel. 

M s is the magnetic energy associated with the cloud. 
Q is the gravitational potential energy. For an axisym- 
metric ellipsoid with concentric density profiles 
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where j3 = arcsine/e is the geometry factor determined 
by eccentricity e = y/l — f 2 and a is determined by the 
density profile as being (1 — a/3) /(I — 2a/5) for a power 
law p cx r~ a and 1.3 1 for an isothermal cloud in equilib- 
rium (jBonnor II1956T) . The intrinsic axis ratio / is smaller 
than the observed axis ratio fobs due to projection. My- 
ers et al. (1991) show that dark cloud cores are more 
likely to be prolate. Condensations in Orion are also sta- 
tistically more consistent with being prolate rather than 
oblate (Li 2002). Based on the more general inversion 
from Fall & Frenk (1983), we derive an averaged true 
axis ratio / corresponding to a sample of a single value, 
namely, the observed f b s for a prolate ellipsoid 

/ = -/ obs ^i(0.5,0.5,-0.5,1.5,l,l-/ o 2 bs ) , (7) 

where T\ is the Appell hypergeometric function of the 
first kind. It is easy to verify that / = 1 for f Q b s = 1 and 
/ = for fobs = 0._ For values in between, we use the 
numerical value of / as a replacement of / in calculating 
the gravitational potential. 

By assuming M. s — and £q =0, a simple virial mass 
can be derived using the virial theorem (Eq. [4]) 
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ap G 
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The virial mass is the minimal mass for a cloud to be self- 
gravitating in the absence of pressure confinement and 



steady magnetic field. A "virial mass ratio" independent 
of B and £ q is defined as 



M 



(9) 



We find that 28 of 30 cores (Fig. [3] and Table [TJ, i.e., 
93% have R mr > 1, 15 cores (50%) fulfil R mr > 5. The 
majority of the core sample are bound by gravity. 

4.2. Stability and Critical Mass 

Now we consider the stability of cores with external 
pressure and the internal magnetic field. The critical 
mass Mq is defined as the maximum mass which can 
be stably supported by internal velocity dispersion and 
magnetic pressure. The two effects can be considered 
separately and then combined in a simple approximation 



M c = Mj + M$ 



(10) 



which is accurate to within 5% co mpared to the more 
rigorous calculations (|McKee 111989ft . 

The Jeans mass for a non-magnetic i sothermal cloud 
(|Bonnor 111951 IMcKee fc Zweibel I fl99l is 
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where c w parametrizes the pressure from turbulence out- 
side of the cloud and c pr parametrize the added internal 
pressure from turbulence's disturbing B field. If we con- 
sider tracers of more diffuse gas in this region, such as 
13 CO, the lincwidth is around 1.5 kms -1 (Melnick et al. 
2011), definitely supersonic. We choose the combined 
factor to be in the middle of the possible range for tur- 
bulent gas cj 2 / c y r 2 ~ 2.1 (jMcKee fc Zweibel l[j"992h . 
The external pressure term can be estimated as 



Pic — Tli 



(12) 



where the mean molecular weight /i is 2.88 and ttih is the 
proton mass. We take a lower bound value of 1 kms -1 
FWHM for estimating Ui C . The inter-cloud density rii C is 
not well known. Radiative transfer analysis of CO, 13 CO, 
and 18o in this region (Li 2002) show that the densities 
are a few times 10 4 cm -3 in the relatively diffuse area. 
We use a lower bound of 10 4 cm~ 3 to calculate the Jeans 
mass given in Tabic [T] The Jeans mass thus derived will 
be a conservative upper limit. 

The maximum mass which can be supported by a 
steady B field alone is 



nBr 2 
G 1 / 2 



(13) 



where c$ ~ 0.12 is given by numerical simulations 
(Mouschovias & Spitzer; Tomisaka, Ikeuchi & Nakamura 
1988). 

Crutcher et al. (1999) detect the Zeeman effect in the 
CN 3mm line near Orion BN/KL. The field strength is 
derived to be 0.19 mG or 0.36 mG, based on different fit- 
ting schemes. This is larger than the B ~ 0.03 mG mea- 
sured in dark clouds. Given the proximity of BN/KL to 
active star formation, the large B could be explained by 
rapid collapses' freezing magnetic flux into high density 
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regions along the line of sight. We thus take a smaller, 
nominal B — 0.1 mG in our calculations of M$. 

Using the critical mass calculated through Eq. [TOl ITTT 
and[T3l we define the "critical mass ratio" 



M 
M r 



(14) 



We find that 26 out of 30 (87%) cores have R c > 1, 22 
cores (73%) fulfil R c > 2. The majority of the Orion 
massive cores are supercritical. 
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Fig. 3. — The mass and size of cores plotted in the same fashion 
as in paper II. The colors of the the open symbols denote the virial 
mass ratio R v i r , with green for a value smaller than 1, red for 
between 1 and 5 and black for > 5. Open stellar symbols are for 
cores containing a protostar and open circles are for starless cores. 
The colors of the solid squares denote the critical mass ratio, with 
green for < 1 and black for > 1. The solid curve is for mass 
supported a steady magnetic field of 100 fiG. The dashed curv e is 
for the empirical threshold for massive star formation (Eg. |15t . 



4.3. Discussion: Is Orion Special? 

Our data suggest mass ratios R V i r > 5 and R c > 1 
for most cores. This is different from other clouds in 
solar neighborhood within ~ 500 pc. In Perseus, for 
example, Foster et al. (2008) and Kirk et al. (2007) find 
an equivalent R V i r < 3. For Ophiuchus, the observations 
of Andre et al. (2007) imply R V i r < 1 based on their 
Table 4 after adopting the dust opacity used here. For 
SCUBAidentified Gould Belt cores, Sadavoy et al. (2010) 
find that only 20 cores out of a total of 354 have M/Mj 
I 2 (using our dust opacities, and excluding their Orion 
data to reduce sample overlap)." Their Mj represents a 
purely thermal Jeans mass, which means that these cores 
will have an even smaller mass ratio using definitions in 
this paper. 

The observed R V i r in Orion may be similar to those 
found in high mass star forming regions, e.g., the two 



clouds studied by Pillai et al. (2011). Although direct 
comparison is hard due to the much larger distance of 
other high star forming regions. 

The calculations differ slightly between studies, e.g., 
in the factors a, f3, and dust opacities. These differences 
alone cannot explain the larger mass ratios observed in 
Orion. Thus, our result that 50% of Orion cores have 
Rvir > 5 and 73% are supercritical is significant. 

The NH3 lines are not well resolved in our data with a 
channel width of 0.62 kms -1 . Using the channel width 
as the upper limit of the intrinsic line width, the veloc- 
ity dispersion of Orion cores cannot significantly exceed 
those of the cores in other regions studied by Kirk et 
al. (2007), Andre et al. (2007), and Foster et al. (2008). 
The high mass ratios is thus a result of higher masses for 
Orion cores. 

Kauffmann et al. (2010) propose an empirical thresh- 
hold for massive star formation 



m lim (r) = 870 M (r/pc) 1 



(15) 



We indeed find 12 cores (40%) above this threshold. The 
excess is not dramatic, but consistent with Orion cores 
being different from other nearby star forming regions in 
them being more massive and more likely to form massive 
stars. 

5. CONCLUSIONS 

We have mapped OMC2 and OMC3 in NH 3 (1,1) and 
(2,2) with both VLA and GBT. The combined single dish 
plus interferometric data provide a rare detailed look into 
the thermal and dynamic properties of a collection of 
massive quiescent cores. Our main results are: 

1. We obtain good temperature measurement for 30 
dust cores. The median core temperature is 17 K. 
The typical uncertainty for derived temperature of 
each pixel is about 2 K. 

2. 12 cores are associated with a protostar. The aver- 
age temperatures of the protostellar and the star- 
less cores are similar suggesting that the heating in 
OMC2-3 region is primarily external. 

3. A total of 28 cores (93%) are gravitationally bound 
(R V ir > 1), and 15 cores (50%) achieve R V i r > 
5. Compared to other Gould Belt clouds, a much 
higher fraction of cores are tightly bounded. 

4. 26 out of 30 cores (87%) are more massive than the 
critical mass defined as the combination of Jeans 
mass and mass supported by a steady magnetic 
field of 100 ^G. 

In summary, this sample of Orion cores, identified in dust 
emission with temperature and turbulence measured in 
NH3 inversion lines, are proven to be well bounded by 
gravity and are mostly supercritical. They will evolve 
rapidly, either collapsing to form a star or fragmenting. 



REFERENCES 



Andre, P., Belloche, A., Motte, F., & Peretto, N. 2007, A&A, 472, 
519 

Benson, P. J. & Myers, P. C. 1989, ApJS, 71, 89 
Bonnor, W. B. 1956, MNRAS, 116, 351 



Chini, R., Reipurth, B., Ward-Thompson, D., et al. 1997, ApJ, 
474, L135 

Crutcher, R. M., Troland, T. H., Lazareff, B., Paubert, C, & 
Kazes, I. 1999, ApJ, 514, L121 



6 



TABLE 1 
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